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We analyze the free energy of a static quark anti-quark pair in quenched QCD at short and large 
distances. From this we deduce running couplings, g 2 (r,T), and determine the length scale that 
separates at high temperature the short distance perturbative regime from the large distance non- 
perturbative regime in the QCD plasma phase. Ambiguities in the definition of a coupling beyond 
the perturbative regime are discussed in their relation to phenomenological considerations on heavy 
quark bound states in the quark gluon plasma. Our analysis suggests that it is more appropriate to 
characterize the non-perturbative properties of the QCD plasma phase close to T c in terms remnants 
of the confinement part of the QCD force rather than a strong Coulombic force. 



I. INTRODUCTION 

On quite general grounds it is expected, that funda- 
mental forces between quarks and gluons get modified at 
finite temperature. More precisely, we expect that forces 
between static quarks, i.e. static test charges in a thermal 
medium, change because the gluons, which mediate the 
interaction between the static quarks, also interact with 
the constituents (quarks and gluons) of the thermal bath. 
In particular, above the deconfinement temperature, T c , 
the potential is expected to be exponentially screened at 
large distances (r ^$> 1/T) [l|. In leading order perturba- 
tion theory this happens due to the generation of a chro- 
moelectric (Debye) mass of order gT (with g being the 
gauge coupling). Beyond leading order, however, chro- 
moelectric and chromomagnetic screening effects cannot 
be separated unambiguously. It is this non-perturbative 
large distance physics, which plays a central role in our 
attempts to understand the bulk properties of the QCD 
plasma phase, e.g. the equation of state and the appar- 
ent deviations from ideal gas behavior found in numerical 
calculations [j|. On the other hand it is the short and 
intermediate distance regime which is most important in 
the discussion of signals which are considered today as 
being suitable to gain information on properties of hot 
and dense matter generated experimentally in heavy ion 
collisions. In this paper we will quantify the temperature 
dependence of the length scale which separates these dif- 
ferent regimes and analyze in detail the properties of the 
QCD coupling constant at short and large distances. 

Although a detailed understanding of screening phe- 
nomena at large distances is still missing, it is evident 
that in this regime the temperature is the dominant scale 
and consequently will control the running of the QCD 
coupling, i.e. g ~ g(T) for (rT » 1, T > T c )0. How- 
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ever, at short distances, r ■ max(T, T c ) <C 1, hard pro- 
cesses dominate the physics of the quark gluon plasma 
even at high temperature and it is expected that a scale 
appropriate for this short distance regime will control 
the running of the QCD coupling, i.e. g ~ g(r). The 
interplay between short and large distance length scales 
plays a crucial role for a quantitative understanding of 
hard as well as soft processes in dense matter. It will, 
for instance, determine the range of applicability of per- 
turbative calculations for thermal dilepton rates or the 
production of jets as well as the analysis of processes that 
can lead to thermalization of the dense matter produced 
in heavy ion collisions. Moreover, the short and interme- 
diate distance regime also is most relevant for the discus- 
sion of in-medium modifications of heavy quark bound 
states which are sensitive to thermal modifications of the 
heavy quark potential as well as the role of quasi-particle 
excitations in the quark-gluon plasma. In all these cases 
it is not immediately evident that temperature is the rel- 
evant scale that controls the running of the QCD cou- 
pling at energies currently relevant in heavy ion physics 
where temperatures may be reached which are only mod- 
erately larger than the phase transition temperature T c . 
An analysis of this question becomes of particular in- 
terest in view of the recently suggested scenario for the 
existence of a large number of Coulombic bound states 
in the QCD plasma phase close to T c 0,|j|. 

It is the purpose of this paper, to firmly establish that 
also at finite temperature the QCD coupling indeed runs 
as function of the length scale r and agrees with the zero 
temperature running coupling at sufficiently short dis- 
tances. In fact, we will show that in the entire regime of 
distances for which at zero temperature the heavy quark 
potential is considered to be described well by QCD per- 
turbation theory [1 H (r < 0.1 fm) the QCD cou- 
pling remains unaffected by temperature effects up to 
T ~ 3T C . Furthermore, we will analyze the interplay 
between the two relevant scales controlling the behavior 
of heavy quark free energies, F(r,T), and we will quan- 



tify the distance scale below which processes in the QCD 
plasma phase are still dominated by properties of the 
QCD vacuum and above which screening dominates the 
physics in the plasma phase. 

We will start in section II with a discussion of heavy 
quark free energies and their relation to the QCD cou- 
pling constant and give some details on our simulation 
parameters in section III. In section IV we discuss the cal- 
culation of running couplings at finite temperature and 
short distances from color singlet free energies and re- 
late these results to properties of the running coupling 
at large distances in section V. Section VI contains our 
conclusions. 



II. FREE ENERGIES AND RUNNING 
COUPLINGS 

A. Free energy of a static quark anti-quark pair in 
lattice QCD 

Our main concern in this study is the determination of 
a running coupling from correlation functions of Polyakov 
loops. In particular we want to make contact to calcula- 
tions performed at zero temperature. Here the running 
of the QCD coupling has been determined successfully 
in lattice calculations from the short distance properties 
of the heavy quark potential 0, 111 UIj an d contact could 
be made withperturbative results obtained in 2-loop cal- 
culations [a, |(| . We thus want to analyze a finite tem- 
perature observable which naturally is related to the zero 
temperature heavy quark potential and has a well defined 
(perturbative) interpretation that allows us to make con- 
tact with perturbative definitions of a running coupling. 
For this reason we will analyze properties of heavy quark 
free energies in the singlet channel. 

Correlation functions of Polyakov loops define the free 
energy of a heavy quark anti-quark pair. One generally 
considers the so-called color averaged free energy, 

e -F(r,T)/T+C = I( TrjL ( x )TrLt( )) , (1) 

9 

where the Polyakov loop, L(x), is defined on lattices with 
temporal extent N T in terms of temporal link variables 
U (x,r) e SU(3), 



where 



L(x) 






(2) 



Furthermore, r = |x| and C is a suitably chosen renor- 
malization constant, which can be determined from a 
matching of finite temperature free energies to the zero 
temperature heavy quark potential [ll], [I3 • 

The color averaged free energy can be considered as 
a thermal average over contributions corresponding to 
quark anti-quark sources in color singlet and octet states, 
respectively, [l3l.ll4j 

1 



-F^r^/T+C 



-Tr(Z(x)£t( )) 



e - Fs {r,T)/T+c = i(TrL(x)TrLt( )} 
8 

iTr(L(x)£t( )) 



(4) 



(5) 



While the color averaged free energy is defined in terms 
of a gauge invariant Polyakov loop correlation function, 
the singlet and octet correlation functions are given in 
terms of a gauge dependent correlator, TrL(x)L T (0), and 
thus have to be evaluated in a fixed gauge. Nonetheless, 
TrL(x)L^(0) is related to an appropriately chosen gauge 
invariant correlator and thus has a proper gauge invari- 
ant interpretation. In fact, when evaluated in Coulomb 
gauge the singlet correlation function constructed from 
the Polyakov loops defined in Eq. [21 may be viewed as 
resulting from gauge fixing non-local but gauge invariant 
operators, i.e. dressed Polyakov loops, where the static 
quark and anti-quark sources are surrounded by gluon 
clouds [Hill. 

For our purpose of defining a running coupling at finite 
temperature the singlet free energy is most appropriate 
as it has at short (r ■ max(T, T c ) <C 1) as well as large 
(rT » 1,T >■ T c ) distances and temperatures a simple 
asymptotic behavior which is dominated by one gluon 
exchange J2J, i.e. 



Fi(r,T) = 



3» 

3-7TT 



9 2 (T) q (T)rT 



3-Trr 



r ■ max(T, T c ) < 1 
(rT»l,T>T c ) 



(6) 



-F(r,T)/T 



-F 1 (r,T)/T 



-F s (r,T)/T 



(3) 



Here we have already anticipated the running of the cou- 
plings with the dominant scales in both limiting regimes, 
although their running, of course, only arises in higher or- 
der perturbative calculations. We also have suppressed 
any additative constants which result from the renormal- 
ization of the free energy [ll| and, in particular, at high 
temperature will dominate the free energy in the large 
distance limit. 

While the relation between the singlet free energy and 
a running coupling is straightforward the definition of a 
running coupling with the help of the color averaged free 
energy is problematic. The exact cancellation of leading 
order perturbative terms in high temperature perturba- 
tion theory, F(r,T)/T ~ (Fi(r, T)/T) 2 , which generi- 
cally does not occur at finite distances and temperatures, 
makes it difficult to define a running coupling which eas- 
ily could be motivated by perturbation theory. While the 
above quadratic relation between F(r,T) and Fi(r,T) 
holds at large distances it has been demonstrated by us 
recently that even at high temperature the short distance 
part of the color averaged free energy is dominated by the 
singlet contribution and F(r, T) ~ Fi(r, T) holds at short 
distances |ll|. In fact, for this reason the determination 
of a screening mass from the exponential fall of F(r, T) 
at large distances, i.e. F(r,T) ~ exp{— fi(T)r}, turned 



out to be quite difficult and strongly dependent on the 
asymptotic form used in fits of the large distance behav- 
ior of F(r,T) 17]. The sub-leading power-like correc- 
tions were found to be strongly temperature dependent 
and turned out to be difficult to control. 



B. Running couplings 

The perturbative short and large distance relations for 
the singlet free energy will be used to define a running 
coupling at finite temperature. In general, the defini- 
tion of a running coupling in QCD is not unique beyond 
the validity range of 2-loop perturbation theory; aside 
from the scheme dependence of higher order coefficients 
in the QCD /3-functions it will strongly depend on non- 
perturbative contributions to the observable used for its 
definition. This is quite apparent when defining the cou- 
pling in QCD cither in terms of the free energy (T = 0: 
potential) 



ay(r,T) = ~F 1 (r,T) 

or its derivative (T = 0: force) 



i(r,T) 



3r 2 dFi(r,T) 
4 dr 



(7) 



(8) 



At low temperature the former necessarily has to change 
sign at some intermediate distance due to the dominance 
of the linearly rising confinement part in the potential 
|18j . The latter, however, stays positive as Fi(r,T) in- 
creases monotonically with r. In fact, for a linear confin- 
ing potential ay defined through Eq.0will become neg- 
ative and drop quadratically while a qq defined through 
Eq. |S|will rise quadratically at large distances. The lat- 
ter gives the possibility of smoothly matching the in- 
creasing coupling in the perturbative regime to the non- 
perturbative increase. This and the poor convergence of 
the perturbative expansion for ay(r) have been reasons 
for analyzing in lattice calculations running couplings de- 
fined through Eq. |H1 (gq-scheme) rather than Eq. (V- 
scheme). We will in the following consider both defini- 
tions as this will also help to distinguish the short and 
large distance regimes at finite temperature. 

In the perturbative regime different definitions of the 
running coupling are uniquely related through the QCD 
/?- function, 



a qq {r,T) = ay(r,T) -r 



day(r,T) 
dr 



ay(r,T) + 9 -^p-f3(g) 



(9) 



with j3(g) = —bog 3 — big 5 — b2g 7 + 0(g 9 ) and universal co- 
efficients b = 11/16tt 2 , & a = 102/(16tt 2 ) 2 . In higher or- 
ders the coefficients of the /3-function are scheme depen- 
dent. At zero temperature the heavy quark potential has 
been calculated in 2-loop perturbation theory [a, Ll • 



From these calculations the 3-loop coefficient, 62, in the 
F-scheme |6j and the gq-scheme Q could be extracted 
and allowed for a detailed comparison of running cou- 
plings determined in perturbative and non-perturbativc 
lattice calculations (Jj, [lfj. Good agreement has been 
found at distances r < O.lfm, which also has been es- 
timated to be the range of validity of the perturbative 
calculations. At larger distances, however, any perturba- 
tively motivated definition of the running coupling will 
also become sensitive to non-perturbative effects and may 
lead to quite different results. 

We will extent here the zero temperature studies of 
the heavy quark potential and the force between static 
charges to finite temperature. In this case the appro- 
priate observable is the heavy quark free energy and its 
derivative. At short distances we follow the approach 
used also at T = and introduce a running coupling by 
analyzing the r-dependence of the force between static 
quark anti-quark sources, Eq. [5J and will compare it 
with the definition of a coupling in terms of the poten- 
tial, Eq. [7| At large distances we will determine a In- 
dependent running coupling directly from fits of F\ (r, T) 
which are motivated by the perturbative large distance 
form given in Eq. We will be more specific on this in 
Section V. 



III. SIMULATION PARAMETERS 

We will analyze in the following properties of heavy 
quark anti-quark pairs in a thermal heat bath of glu- 
ons, i.e. we consider correlation functions of Polyakov 
loops in the SU(3) gauge theory (quenched QCD) at fi- 
nite temperature calculated on Euclidean lattice of size 
N 3 x N T . All our calculations have been performed on 
lattices with spatial extent N a = 32 and N T — 4, 8 and 
16 using the tree level Symanzik-improved gauge action 
(19J . It has been verified by us earlier [llj that this choice 
of action and lattice parameters is sufficient to suppress 
finite volume effects and finite size effects such as the 
breaking of rotational symmetry in the analysis of cor- 
relation functions at short and intermediate distances. 
Our simulations have been performed in the tempera- 
ture range T c < T < 12 T c . The simulation parameters 
are summarized in Table [I] The temperature scale for 
the Symanzik-improved action was obtained earlier from 
calculations of the string tension at zero temperature and 
a determination of the critical coupling for the deconfine- 
ment transition on lattices with temporal extent N T = 4 
and 6 ;20]. 

The calculation of singlet free energies has been per- 
formed in Coulomb gauge, which on the lattice is realized 
3 

by maximizing Tr > U^x., r) in each time slice, i.e. for 

fixed t. A residual gauge degree of freedom is fixed by de- 
manding N^ Uo{x,t) to be independent of r. Typically 





/3 


T/T c 


N r = 4 N T = 8 N T = 16 


1.03(1) 


4.090 


- 


1.05(2) 


4.100 


4.5592 


1.10(1) 


4.127 


- 


1.15(1) 


4.154 


- 


1.20(3) 


4.179 


4.6605 


1.24(1) 


4.200 


- 


1.30(1) 


4.229 


4.7246 


1.50(3) 


4.321 


4.8393 5.4261 


1.60(2) 


4.365 


4.8921 


1.68(2) 


4.400 


- 


2.21(5) 


4.600 


- 


3.0(1) 


4.839 


5.4261 


6.0(3) 


- 


6.0434 


9.0(3) 


- 


6.3910 


12.0(5) 


- 


6.6450 



TABLE I: Parameters of the simulations on 32' x N T lattices 
using the tree-level Symanzik-improved action. 



we have analyzed the correlation functions on 100-500 
independent gauge field configurations. 

As we are interested in the short distance behavior of 
the heavy quark free energy it is important to correct for 
the violation of rotational symmetry which is most pro- 
nounced in this region. Following [lOj we have replaced 
F(r,T) by F(ri,T) (similarly for F\ and F 8 ) where 77 
relates the separation between the static quark and anti- 
quark sources to the Fourier transform of the tree-level 



lattice gluon propagator, D 



I IV: 



= 47T 



w d 3 k 
« (2tt) 3 



exp(ifc • r)Doo(fe) 



(10) 



For the Symanzik-improved action the time-like compo- 
nent of Dfn, is given by 



D^{k)=A 



3 

E 

i=l 



■ 2 ^i 

sin — 
2 



1 • 4 
-sm 



(11) 



This procedure removes most of the lattice artifacts as is 
evident from the smooth short distance behavior of the 
free energies shown in Fig. ^ Therefore in what follows 
we will always show free energies plotted versus 77, but 
will suppress the subscript I. 



IV. RUNNING COUPLING CONSTANT AT 
SHORT DISTANCES 

In Fig. n we show our results for the singlet free ener- 
gies for several temperatures above T c . The free energies 
have been calculated in Coulomb gauge and have been 
renormalized by matching the short distance part to the 
zero temperature heavy quark potential of Ref . |l(J . 
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FIG. 1: Heavy quark free energies in the singlet channel cal- 
culated for several values of the temperature above T c (a) and 
the short distance part of the singlet free energy at T/T c — 1.5 
(b). The dashed line shows V s tring(r) defined in Eg. 1131 and 
the solid line is the perturbative result of Ref. 10] for V q q(r). 
This has been matched smoothly to V 3 t r ing(r) at r ~ 0.1 fm. 
Physical units have been obtained by using ^J~o = 420 MeV. 



Fig. ^ shows the short distance part of the singlet 
free energy calculated at T = 1.5T C on lattices with 
temporal extent N T = 4, 8 and 16. This corresponds 
to lattice spacings ranging from a ~ 0.12 fm down to 
a ~ 0.03 fm. Apparently the short distance part of the 
singlet free energy agrees quite well with the zero temper- 
ature heavy quark potential including a perturbatively 
calculated Coulomb term Hi, 



V m (r) 



Aay(r) 
"3 r 



(12) 



and shows no significant cut-off dependence. On the 
other hand, deviations from a confinement potential with 
a constant Coulomb like term that arises from universal 
string fluctuations, 



V t 



string 
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err 
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FIG. 2: The running coupling in the gg-scheme determined 
on lattices of size 32 3 x N T with A'',- = 4 (open symbols) and 
8 (filled symbols) from derivatives the short distance part of 
the singlet free energy (T = 0: from the force) at different 
temperatures. The relation of different symbols to the values 
of the temperature are as in Fig. la. The various lines are 
explained in the text. 



FIG. 3: The running coupling in the V-scheme determined 
on lattices of size 32^ x N T with N T — 4 (open symbols) 
and 8 (filled symbols) from the short distance part of the 
singlet free energy (T = 0: from the potential) at different 
temperatures. The relation of different symbols to the values 
of the temperature are as in Fig. la. The dashed (solid) lines 
for Q qq (ov) represent the zero temperature results of Ref. |l0|) 
which are continued to shorter distances using our new results. 



are clearly visible at these short distances. This already 
indicates that the short distance behavior of the singlet 
free energy is consistent with a running coupling that 
is controlled by the quark anti-quark separation, r, and 
shows no or only little temperature dependence. We 
also note that after having renormalized the free energy 
through a matching at short distances the large distance 
behavior is completely fixed. Above T c the singlet free en- 
ergy approaches a temperature dependent constant value 
at large distances which changes sign for T ~ 3 T c . 

To analyze the T and r-dependence of the coupling we 
first follow the approach used at T = |3, Q and define 
a qq (r,T) through Eq. |SJ The derivatives of the singlet 
free energy with respect to the distance, di*i(r, T)/dr, 
are obtained from a finite difference approximation using 
results at neighboring distances. We compare our finite 
temperature results to the high statistics calculation per- 
formed at zero temperature [lfj in Fig. [21 In this figure 
the results obtained from the numerical calculation of the 
heavy quark potential at T = and distances r > 0.1 fm 
are summarized by a fat black line. At shorter distances 
a thin line represents the result of a perturbative calcu- 
lation of the force [2, Il3 which is based on the 2-loop 
calculation of the heavy quark potential |a, Q . This per- 
turbative result is smoothly matched to the lattice data 
at r ~ 0.1 fm. Also shown in the figure as a dashed line 
is the effective coupling extracted from the confinement 
potential V s tring using Eq. |SJ It agrees quite well with 
the lattice data for r > 0.1 fm but shows strong devia- 
tions from the perturbative as well as lattice calculation 
at shorter distances. 

Our numerical results on a qq at distances smaller than 
0.1 fm cover also distances substantially smaller than 



those analyzed so far at T = 0. They clearly show the 
running of the coupling with the dominant length scale 
r also in the QCD plasma phase. For temperatures be- 
low 3T C we find that a qq agrees with the zero temperature 
perturbative result in its entire regime of validity, i. e. for 
r < 0.1 fm. At these temperatures thermal effects only 
become visible at larger distances and lead, as expected, 
to a decrease of the coupling relative to its zero tem- 
perature value; for larger temperatures thermal effects 
influence also the short distance behavior at distances 
r < 0.1 fm. 

At distances larger than r ~ 0.1 fm non-perturbative 
effects clearly dominate the properties of a qq . It thus is 
to be expected that the properties of a running coupling 
will strongly depend on the physical observable used to 
define it. To quantify this we analyze directly the short 
distance behavior of the renormalized singlet free ener- 
gies and define a v (r,T) through Eq. At T = the 
singlet free energy simply is the heavy quark potential 
which for distances r > 0.1 fm is quite well described by 
V s tring (r). In the case of the string potential, V s t r ing, one 
would obtain with the perturbative ansatz given in Eq. 
for the running coupling ay(r, 0) = 7r/16— 0.75ot 2 , which 
changes sign at r ~ 0.25 fm. We expect to find a similar 
behavior also when using F± (r, T) at temperatures close 
to T c . As can be seen from Fig. such a behavior is 
indeed found for T c < T < 3T C . This again reflects the 
importance of remnants of the confining force in the QCD 
plasma phase. For larger temperatures ay stays positive 
reflecting the fact that F\ (r, T) approaches a negative 
constant at large distances. We note that at short dis- 
tance ay(r,T) > a qq (r,T) as expected from the pertur- 



bative relation, Eq. between both definitions of the 
running coupling. At distances r > 0.1 fm the couplings, 
however, merely reflect the non-perturbative properties 
of the observable used to define them. 

Let us return to an analysis of the properties of the 
running coupling defined in the gg-scheme. At temper- 
atures close to T c the coupling ay(r, T) stays close to 
the zero temperature value up to distances r ~ 0.3 fm. 
At these distances the strong r-dependence of a qq (r,T) 
mimics the linear rising part of the zero temperature 
confinement potential. Remnants of the confining force 
thus survive the deconfinement transition and play an 
important role in quark anti-quark interactions at in- 
termediate distances up to r ~ 0.3 fm. At distances 
larger than ~ 0.5 fm these are, however, rapidly screened 
also at temperatures close to T c . We use the maxima 
in a qq (r,T) to define a length scale r scrcon which sepa- 
rates the short distance regime from the large distance 
regime. This is shown in Fig. 0] In the entire tem- 
perature interval analyzed by us we find that r scrccn is 
inversely proportional to the temperature, i.e. we find 
Tscrcon = (0.48 ± 0.01)fm • T c /T. 

The fact that a qq (r,T) can become large at some dis- 
tance also in the deconfined phase of QCD has recently 
been exploited to discuss the scenario of a strongly in- 
teracting fluid of quasi-particles describing the thermo- 
dynamics above but close to T c 0, 3 • Our analysis of 
the running coupling shows that up to a certain distance 
scale confining features of the heavy quark potential in- 
deed survive in the plasma phase and thus may support 
such a scenario. We note, however, that this effect is not 
related to an unexpectedly large coupling arising from 
thermal effects but on the contrary to the survival of vac- 
uum physics below a certain characteristic length scale. 
In particular, there is no evidence for a larger coupling 
in the Coulomb part of the heavy quark potential arising 
from thermal effects. From this point of view it may be 
questionable whether the plasma phase really can sup- 
port the existence of well localized colored bound states 
as it is advocated in [3|. 



V. T DEPENDENCE OF THE COUPLING AT 
LARGE DISTANCES 

We now turn to an analysis of the large distance struc- 
ture of heavy quark free energies and discuss the deter- 
mination of a temperature dependent coupling from it. 
As inferred from the r-dependence of the coupling the 
crossover from the short to large distance regime sets 
in rather abruptly. In particular, the screening of the 
heavy quark free energies leads to an exponential sup- 
pression of a qq . In order to extract the T-dcpendcncc 
of the QCD coupling conventionally used to describe the 
large distance properties of QCD at high temperature 
this screening effect should be eliminated. In Fig. we 
show r(F'i(oo, T)—Fi(r, T)) on a logarithmic scale. Aside 
from deviations at short and intermediate distances this 




FIG. 4: The location of the maximum in a qcl (r,T) at fixed 
temperature versus temperature in units of the deconfinement 
temperature. Open (filled) symbols correspond to calcula- 
tions on lattices with temporal extent iV T = 4 and 8, respec- 
tively. The solid line corresponds to r scrccn = 0.48fm • T c /T. 




FIG. 5: The singlet free energy versus rT for several temper- 
atures above T c obtained from calculations on lattices with 
temporal extent N T — 4 (open symbols) and N T = 8 (filled 
symbols) . 



is seen to decay exponentially. We fit the large distance 
part of F\ (r, T) with an ansatz motivated by the Debye 
screened perturbative result, 



Fm (r,T) 4a(T) 



T 



3rT 



exp{- ^4ma{T)rT} + c(T) , (14) 



where a(T) and a(T) are used as two independent fit 
parameters, which at large temperature, i.e. in the per- 
turbative limit, are expected to coincide, 



lim a(T)/a(T) 

1 — >oc 



1 



(15) 



We note that the coupling a(T) introduced here is closely 
related to the definition of ay given in Eq. {7\ 



Both fit parameters as well as their ratio are shown in 
Fig. [fj] Within statistical errors the above limit indeed 
is reached for temperatures T/T c ~ 6. At smaller tem- 
peratures, however, deviations from unity are large and 
reach a(T)/a(T) ~ 5 close to T c . Similar to the ambigu- 
ities that exist for the definition of a running coupling at 
large distances this suggests that for not too large tem- 
peratures, T c < T < 2T C , any definition of a tempera- 
ture dependent running coupling will strongly depend on 
the physical process used to define a(T) or equivalently 
g(T). Nonetheless for T > QT C the temperature depen- 
dence seems to follow the logarithmic behavior expected 
from perturbation theory. This is shown in the figure by 
the solid lines with the dotted error band which represent 
an appropriately rescaled perturbative running coupling, 
a(T) — 2.095(82)a pcr t(T) where we have used for a pcr t 
the 2-loop perturbative running coupling, 
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and relate A™- to the critical temperature for the de- 
confinement transition, T c /Aj^g ~ 1.14(4) H,0|. The 
rescaling factor has been determined from a common fit 
of the data for a(T) and a(T) at temperatures T > 6 T c . 
At least in this temperature regime all higher order per- 
turbative as well as non-perturbative effects seem to be 
well described by a rescaling of the coupling with con- 
stant factor. A similar observation has been made for 
screening masses determined in an SU(2) gauge theory 
at high temperature |22|. 

Also shown in the upper part of Fig. is the max- 
imal value of the running coupling determined from 
the short distance behavior of the singlet free energy, 
<2qq(r SC rccn, T). As can be seen, this coupling is signifi- 
cantly smaller than a(T) determined as coupling strength 
of the Debye screened Coulomb potential which describes 
the long distance part of the singlet free energy. We stress 
again that a(T) found here for T < 3T C at large distances 
is not appropriate to characterize the Coulombic part of 
F\ (r, T) at short distances. As discussed in the previous 
section this is still controlled by an almost temperature 
independent coupling ay (r, T) . 



VI. SUMMARY 
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We have performed a detailed study of the singlet free 
energy of a static gq-pair in quenched QCD (SU(3) gauge 
theory) at short and large distances. We have shown 
that at sufficiently short distances the free energy agrees 
well with the zero temperature heavy quark potential 
and thus also leads to a temperature independent run- 
ning coupling. This short distance regime is temperature 
dependent and reduces from r ~ 0.5 fm at T ~ T c to 
r ~ 0.03 fm at T ~ 12 T c . At high enough tempera- 
tures, T > 6T C , the large distance behavior of the free 
energy is qualitatively similar to what is expected in per- 



FIG. 6: The temperature dependent running coupling de- 
termined from the large distance behavior of the singlet 
free energy on lattices with temporal extent N T — 4 (open 
symbols) and N T — 8 (filled symbols). The upper figure 
shows a(T) = g (T)/Atv (dots) and the value a qq (r SC iccn,T) 
(squares) determined from the short distance behavior of the 
singlet free energy (see Fig. 3). The figure in the middle shows 
5(T) = g 2 (T)/4iv and characterizes the temperature depen- 
dence of the screening mass. The lower figure gives the ratio 
of both fit parameters. The solid lines with the dotted error 
band are discussed in the text. 



turbation theory and allows for a consistent definition of 
a temperature dependent running coupling characteriz- 
ing large distance properties of QCD thermodynamics. 
However, at temperatures close to T c the definition of a 
temperature dependent running coupling is not unique. 
Our analysis suggests that it is more appropriate to 
characterize the non-perturbative properties of the QCD 
plasma phase close to T c in terms remnants of the con- 
finement part of the QCD force rather than a strongly 
coupled Coulombic force. 
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